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Abstract. The influence of time dependent medium modifications of low mass 
vector mesons on dilepton production is investigated using nonequilibrium 
quantum field theory. By working in the two-time representation memory 
efi^ects are fully taken into account. For different scenarios the resulting dilep- 
ton yields are compared to quasi-equilibrium calculations and remarkable dif- 
ferences are found leading to the conclusion that memory effects can not be 
neglected when calculating dilepton yields from heavy ion collisions. Emphasis 
is put on a dropping mass scenario since it was recently claimed to be ruled 
out by the NA60 data. Here the memory effects turn out to be particularly 
important. 
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1. Introduction 



Relativistic heavy ion reactions allow for studying strongly interacting matter at 
high densities and temperatures. One of the main objectives is the creation and 
identification of new states of matter, most notably the quark-gluon plasma (QGP). 
Photons and dileptons do not undergo strong interactions and thus carry undis- 
torted information on the early phases of the fireball because the production rates 
increase rapidly with temperature. Photon spectra are a suitable observable for 
the temperature whereas in the low mass region dileptons couple directly to the 
light vector mesons and reflect their mass distribution. They are thus considered 
the prime observable in studying mass (de-)generation related to restoration of the 
spontaneously broken chiral symmetry. In order to draw conclusions from the data 
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(eg. CERES [ 1 , NA60 [H E] or HADES) a precise theoretical description of the 
medium effects has to be estabhshed. 

During the early stages of a heavy ion reaction the system is out of equilibrium. Also 
at later times, while the system further expands, it is (slightly) out of equilibrium. 
In standard settings the system is thought to expand through quasi-equilibrium 
within a fireball scheme: The dilcpton production is then calculated under the as- 
sumption of instantaneous forms of the production rate at the given temperatures 
and convolution with the expansion ['51, 151. 

In contrast, we account for possible retardation and off equilibrium by using a 
nonequilibrium quantum field theoretical description based on the Kadanoff-Baym 
equations [IHI. We simulate modifications of the light vector mesons by introduction 
of a certain time dependence of the self energies, which then allows for analyzing the 
dynamics of the mesons' spectral properties as well as that of the resulting dilepton 
rates. In [ |7] we introduced time scales which characterize the memory of the 
spectral function and compared them to typical time scales in heavy ion reactions. 
We found that changes are not generally adiabatic and memory effects have to 
be included. Here we present dilepton yields calculated for constant temperature 
to give a first impression of the memory effects. Furthermore we concentrate on 
dilepton yields for a dropping mass scenario within a heavy ion collision like that 
measured by the NA60 experiment Comparison of the experimental data with 
equilibrium calculations has lead many people to conclude that Brown- Rho scaling 
is ruled out by the data (see [El El for a contrary opinion). We show that for this 
case memory effects are especially important and before drawing any conclusion 
they should be incorporated into any calculation. 



2. Equilibrium vs. Nonequilibrium 

We compare the usual equilibrium approach, in that the dilepton rate is given by 
the well known formula 
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nonequilibrium formalism, in which the propagators of the p-meson and the virtual 
photon are calculated using the general nonequilibrium formulas and the dilepton 
rate for a time dependent but otherwise homogeneous situation is given by [|Z|: 
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with the transverse virtual photon propagator D^j,, following the generalized fluc- 
tuation dissipation theorem 



ret < T^adv 



(3) 
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with all integrals, time variables and further indices implicit. In this case the 
£)ret/adv ^]^g £j,gg retarded and advanced virtual photon propagators while = 
n< , the photon self energy, which follows via VMD by 



Note that the couplings are to be taken at the correct vertices, separated in time. 
The meson propagator -D^^^ also follows (|2J), with g< = the meson self energy, 
and the dressed retarded and advanced meson propagators which follow the equation 
of motion that in a spatially homogeneous and isotropic medium is given by 

{-dl ~mf~ k^) i^-t (k, hM) - r dtY^^TiK h,l)Dl%{]^, t, h) - 5{t, - h), 

Jt2 

(4) 

and Z)^'^(k, ti,t2) = Dp°-^(k, t2,ti). The integral over past times in Eq. (0) as well 
as the time integrations in Eq. have encoded the finite memory of the system. 
In [[7] we show that for describing a heavy ion reaction it is essential to retain this 
dynamical information, because timescales for the adaption of the meson's spectral 
properties to the evolving medium are of the same order as the lifetime of the re- 
garded hadronic system. Hence, an assumed adiabatic instantaneous approximation 
as in Eq. (QJ can not properly describe the situation in a heavy ion collision, and 
remarkable differences in the resulting calculated yields occur. 



3. Nonequilibrium dilepton production from evolving media 

3.1. Self energies 

The dynamic medium evolution is introduced by hand via a specified time dependent 
mass and retarded meson self energy I1'^°^{t,uj) [|7]. From that the self energy S< 
follows by introducing a background temperature of the fireball. The fireball is 
assumed to generate the time dependent self energy S''*'* and, assuming a quasi 
thermalized system, the p-meson current-current correlator is given by 

E<T(r,u;,k) = 2inB(T(r))ImS-*(r,c^,k), 

which follows from the KMS relation, being valid for thermal systems. The latter 
is a rather strong assumption, but necessary in order to proceed. 

The medium effects are introduced via a specific evolving retarded self energy 
of the vector meson. A simple self energy 

ImS'-<='(tj,r) = -cjr(T), (5) 

with a k- and oj-independent width F, leads to a Breit-Wigner distribution for the 
spectral function. The time dependence is being accounted for by introduction of 
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the parameter t. For the k = mode, the full self energy for coupling to = | 
-resonances is given by 

J ^rct. 1 ^('^) f [rnpY uj^ETRir) 

ImL (T,ui,k. — 0) — — -] qi „ , ,o ^^^[t)- 

(6) 

with E = \J m?j^ + k'^ ~ tiin and and m^r the masses of the resonance and the 
nucleon respectively. T ^ is the width of the resonance and gi the isospin factor [ 
I10| . Other, more involved microscopic forms of the self energy Iml]''°* can also be 
invoked straight forwardly, see [Ej. 



3.2. Quantum interference 

If the systems properties change fast enough, as in relativistic heavy ion collisions, 
the quantum mechanical nature of the regarded systems can lead to oscillations 
and negative values in the changing spectral functions, occupation numbers and 
production rates as well as interferences that one does not get in semi-classical, adi- 
abatic calculations. The rate calculated here possesses the full quantum mechanical 
information incorporated and contains "memory" interferences that might cause 
cancellations - hence the rate has to be able to become negative while the time 
integrated yield always stays positive as the only observable physical quantity. An 
intriguing example for the occurring oscillations in electron-positron production is 
shown in Figs. ^ and [21 
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Fig. 1. Rate for the change of the 
mass from m = 400 MeV to 770 MeV 
(constant r=150 MeV and constant 
T = 175 MeV) directly after the self 
energy has reached its final form (af- 
ter 7.18 fm/c) and 2 fm/c later. 



Fig. 2. Yield corresponding to rates 
in the left plot integrated over time in- 
tervals starting with starting changes 
and ending after the indicated dura- 
tion. Also yields at earlier and later 
times are shown. 
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3.3. Yields at constant temperature 

In order to get a first impression of how tlie memory effects affect dilepton yields, 
we investigate dielectron- yields from systems witli a constant size at constant tem- 
perature. Therefore we apply linear changes to the system over a duration of 7.18 
fm/c. The results for the yields integrated over this time for four different scenarios 
are shown in Fig. O below. The largest effect is found for the mass shift of the 
/5-meson: The yield is increased by a factor of about 1.8 in the range from 200 
to 450 MeV due to memory effects. The other cases show differences, but not as 
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Fig. 3. Comparison of the dynamically calculated (solid) to the Markovian (dashed) 
dilepton yields from an interval of duration Ar = 7.18 fm/c in that the self energy was 
changed linearly as indicated in the corresponding figure. Initial quantities have index 1, 
final ones index 2. The temperature was kept constant at T = 175 MeV. 

pronounced - we see a minor enhancement of the lower mass tail due to the Bose 
factor in all cases. In the case of a broadened p-meson, one can nicely see that the 
dynamically calculated yield is a broader distribution and that for the case of the 
coupling to the N(1520) without broadening the resonance peak is stronger, due to 
the system's memory of it, and the vacuum peak is shifted further to the right due 
to memory of the level repulsion effect. 

To conclude this section we state that for the p-meson we found exactly the mod- 
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ifications of the yields that one would expect qualitatively when including memory 
effects for the constant temperature case. In a dynamical fireball setup with full de- 
coupling in the end we find similar results but more pronounced differences between 
the quasi-equilibrium and nonequilibrium calculations [[?]. 

4. Dilepton yields from heavy ion reactions 

Triggered by the recent discussion on whether Brown- Rho scaling [ is 'ruled 
out' by the NA60 data, presented at Quark Matter 2005 [E], we concentrate on 
this issue and calculate dilepton yields for both suggested mass modifications of 
the p-meson within the nonequilibrium formulation, which has been introduced in [ 
U]. First quantitative but simplified calculations within this work have shown that 
memory effects, which are neglected in equilibrium calculations (recent calculations 
for the NA60 scenario can be found in [^] and have an important influence 
on the final dilepton yields, especially for dropping mass scenarios. An analysis of 
which scenario describes the data correctly demands the proper inclusion of memory 
effects. The validity of equilibrium calculations will depend on the strength of these 
effects. The two mass parameterizations we use are the temperature and density 
dependent one used by Rapp et al. [ lOj 

m; - mp(l - 0.15nB/7io) [l - {T/T^ff''' , (7) 

assuming a constant gauge coupling g in the vector meson dominance (VMD) - 
coupling, and one motivated by an improved version of Brown- Rho scaling [El El j 
where 

TO* = TOp(l - O.lSrtB/no), (8) 
with a modified gauge coupling g*, in such a way that g* is constant up to normal 
nuclear density tiq, while from then on m*p/g* is taken to be constant [ I14|. This 
is due to the fact that VMD, which accounts for most of the shape of the dilepton 
spectrum, is violated for most temperature and density regions (see [El El and [ 
I15p. Brown and Rho argue [El El that because lattice calculations tell us that 
the pole mass of the vector mesons does not change appreciably up to T = 125 
MeV, the parameterization using [1 — (T/Tc)"]'* (with positive d and integer n) 
overestimates the mass shift. They find temperature dependent effects to be an 
order of magnitude smaller than the density dependent effects, and hence suggest 
to concentrate on the density dependent part. Furthermore they point out that due 
to the violation of VMD the overall dilepton production in dense matter should be 
reduced by a factor of 4 compared to Rapp's calculations, which we do not take 
into account here. Following [ITBl 13, we use a fireball model with a cylindrical 
volume expansion in the ±z direction: 

l/(r) = (zo + «.T)7r(ro -1-0.5 air2)2, (9) 

where zq is equivalent to a formation time tq = 1 fm/c, = c is the primordial 
longitudinal motion, 7'o ~ 5.15 fm is the initial nuclear overlap radius, and a±^ = 



Nonequilibrium dilepton production 



7 



0.08 c^/fm is the radial acceleration. We start the calculation at T = Tc = 175 MeV 
and set the meson self energy to zero for temperatures below chemical freezeout 
at T = 120 MeV. The hadronic phase lives for about 6 fm/c. The calculation is 
done for each momentum mode between fc = and k ^ 1.5 GeV and integrated 
over momenta. 

The results for the dimuon spectra are shown in Figs. 0]and|Hl We find a strong 
enhancement and also modification of the shape in both cases where for the density 
and temperature dependent parameterization the dynamic calculation yields about 
a factor of 4 more dileptons around an invariant mass of 500 MeV as compared to 
the static case, while the other parameterization shows more enhancement x3) 
around 650 MeV. 

The reason for the strong differences are the following: First, effectively the meson 
mass approaches its vacuum value more slowly in the dynamic calculation. Lower 
masses are enhanced due to the Bose factor and if the spectral function has a 
lower mass for a longer time this means more enhancement. Furthermore does the 
memory of higher temperatures increase this enhancement for all masses. Finally, 
the modified coupling in Eq. ^ suppresses the early contributions. That means 
that for the static case all low masses in the spectral function and high temperatures 
are suppressed whereas in the dynamic case these low masses and high temperatures 
still contribute when the couling is larger and hence they are less suppressed by the 
VMD-coupling. 
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Fig. 4. Dimuon yield from decay- 
ing jO-mesons calculated out of equi- 
librium and within the usual instan- 
taneous approximation using param- 
eterization (0) . A difference of about 
a factor of 4 is clearly visible in the 
regime around 500 MeV. 



Fig. 5. Dimuon yield from decay- 
ing p-mesons calculated out of equi- 
librium and within the usual instan- 
taneous approximation using param- 
eterization The most significant 
difference (about a factor of 3) is vis- 
ible between 600 and 700 MeV. 
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5. Summary and Conclusions 

We calculated dilepton yields within a non-equilibrium field theory formalism, based 
on the KadanofF-Baym equations. We investigated possible medium modifications of 
the p-meson for constant temperature and 'Brown-Rho'-scaling within an evolving 
fireball. Special attention was put to possible retardation effects concerning the 
off-shell evolution of the vector mesonie excitations. 

The full quantum field theoretical treatment leads to oscillations in all men- 
tioned quantities when changes in the self energy are performed. This oscillatory 
behavior reveals the quantum mechanical character of the many particle system, 
present in the investigated heavy ion reaction. 

Comparison of dynamically calculated yields with those calculated assuming 
adiabaticity reveals moderate to strong differences. For the dropping mass scenarios 
investigated, an enhancement of about a factor of 3-4 below the p-vacuum mass 
and significant differences in shape have been found, leading to the conclusion that 
memory effects must not be neglected in precision calculations of dilepton yields 
from relativistic heavy ion collisions. 
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